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Abstract 



Oh! 

r-| , We investigate the CP violating form factor of the ZZZ and ZZ^y ver- 

tices in the pair production of Z° bosons. Useful observables in aziniuthal 



distributions are constructed to probe CP nonconservation which may 
originate from these vertices. A simple Two Higgs Model of CP violation 
is used as an illustration. 
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In the near future, with the availabihty of experimental data at energy around the 
electroweak breaking scale, one expects to learn about the structure of the cubic and quartic 
self-interactions of gauge bosons. So far, these interactions have not been directly tested in 
any experiments. 

One exciting possibility is that such interactions will give new insights on CP violation, 
whose physical origin has not been understood with satisfaction yet. The observation of CP 
violation in the kaon system can be explained in various ways within the framework of gauge 
theories, and choosing between them requires additional observation of CP violation. With 
this in mind, it is interesting to look for CP violating signals which may be induced by the 
self-interactions of gauge bosons. We discuss here one such possibility, where the coupling of 
three neutral gauge bosons has a CP violating term in it. We first did a model-independent 
discussion based on the most general form factors. Then a simple model, the two Higgs 
doublet model, is used as an illustration of how the form factors may arise in a realistic CP 
violating theory. 

I. HELICITY AMPLITUDES. 



Such CP-odd term is indeed allowed in general on fundamental grounds, as is obvious 
from the general parametrization of the cubic coupling of gauge bosons [|l|-0]. Most theoret- 
ical studies along this direction have been done l^j-^ only for the process e^e+ -^ W~W^ . 
The effect of CP violation in e~e^ — >■ Z^Z^ has not been thoroughly carried through |0] 
and there is a need of detailed analysis. This motivates us to perform a careful model- 
independent study. In Section IV, a simple Two Higgs Model is used as an illustration. 
We follow the helicity formalism for the Z^ pair production, e~(cr)e"^(o") -^ Z^{X)Z^{X), 
outlined in Appendix D of Ref . ^ . Here we include explicitly effects from the form factors 
fi and /s which describe the vertex V{P) — >• Z{q)Z{q') for out-going on-shell Z° bosons, 
where the incoming particle V is either another Z-boson or a photon: 
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where s = P'^. Note that f^ term is CP-odd. The /s term, although CP-even, is included 
for completeness. The helicity amplitudes are given by 
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The kinematic variables are defined as usual, 7"^ = 1 — /?^ = AM^/s. The amplitude for 
the initial helicity configuration a = a is highly suppressed due to helicity argument in 
the high energy limit ^/s ^ mg. Therefore we are only interested in the cases for which 
a) = ±1. The relevant Wigner d functions appearing in Eq.(||) are listed below: 



Aa=^ia 



^?,±2(0) = -rf-i,T2(0) 



±-(1 ± cos 0) sine. 



<±i(0) = rfii,Ti(0) = o(i±cose), 
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dUQ) = -rf-in(0) 
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sine. 



In the standard electroweak model at the tree level, the elements Ax^x'{Q) come from the 
t-channel exchange diagram. The electron couplings g^a to the Z^ boson are specified by 
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After simplification, we summarize the result for various cases AA = A — A' as follows. 
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When the kinematic variables of the two identical Z^ boson are interchanged, i.e., 

(A,A') ^(A',A), e^TT-e, <l>^7r + $, (6) 

the amplitude is unchanged because of the Bose symmetry, if one includes a negative sign 
coming from the azimuthal $ rotation exp(iAcr7r). 
The usual CP transformation is 



(A,A') ^(-A,-A'), e^TT-e, $^7r + $. 



(7) 



However, we can simplify this CP transformation by incorporating the Bose symmetry in 
Eq.(BD. The resulting CP transformation becomes 



(A, A') -^ (—A', —A), e, $ unchanged. 



(8) 



The situation now becomes very similar to our previous analysis p in the process e e^ — >■ 
W-W+. 

If CP is conserved (when /4's are turned off), we have the following relation for the 
amplitudes in our phase convention: 



-M.,^;A,A'(0)=-^<x,^;-A',-A(e). 



(9) 



This equality will be destroyed by the presence of CP violating form factors /4 in channels 

(A,A') = (0,±)or(±,0). 



II. SPIN-DENSITY MATRICES. 

To avoid studying the complicated event topology in the 4-ferniion final configuration 
from the decays of the Z^ pair, we concentrate our attention to the decay of a single Z°. 
This strategy is equivalent to the study of the density matrix for one of the Z^ bosons. 

We only look at the Z^ boson at the scattering angle and temporarily ignore the 
recoiling one, which is considered as produced at the scattering angle tt — 0. The polar 
angle ip and the azimuthal angle (p are defined in the Z° rest frame for the lepton £~ in the 
decay Z^ — * £^£^. We define the axes of the rest frame of Z° as follows. The z-axis is along 
the direction of motion of Z° in the e~e^ cm. frame. The x-axis lies on the reaction plane 
and toward the direction where G increases. The y-a.xis is given by the right-hand rule. 

The angular distribution of i~ from the Z" — >■ i~i^ decay is specified by the the spin 
density matrix pij of the Z^ boson. 

p(e),, = AfiQ)-' E a^.,.;.a'(0)-m:,.,,a'(0) • (10) 

cr,a-,A' 

Here A/" is the normalization such that Trp = 1. p is hermitian by definition. The normalized 
distribution for i~ is given by 

^^^ = ii ^ ^4(1 + /icos^)V(e)++ + (1 - /icos^)V(e)__ + 2p(e)oosin2^ 

— 2v2Re p(9)+,o(l + hcosijj) sin-?/'cos0 + 2v2Imp(6)+,o(l + hcosip) sin'?/'sin0 

— 2v2Re p(0)_^o(l — hcosip) sin-?/;cos0 — 2v2Imp(0)_^o(l ~ hcosip) sin'?/'sin0 

+2Rep(e)+_(l-cosV)cos20-2Imp(0)+_(l-cosV)sin2(/ 

(11) 

The two contributions come from helicity configurations iJiih = 1) and i^ih = —1), with 
different weights, 

w^=gl/{gl + gl), w+ = gl/{gl + gl) , w.+w+ = l. (12) 

In our present phase convention, if CP were conserved {i.e. when /4 = 0), we would have 
the following identities. 

p(0)a,a' = p(7r - e)_A,-A' , (13) 

based on the transformation in Eq.(|^). Similar expressions were first noticed in Ref. [Q on 
the process e~e^ -^ W~W^ and in Ref. P] on the process e~e^ -^ ti. 



III. CP VIOLATING OBSERVABLES 

Under CP conjugation, we change variables 0— s^tt — O, ■i/'^Tr — ■?/', and -^ —0. 
The distribution in Eq.(|TT]) is transformed into itself if we assume CP conservation as in 



Eq. (|T3]) . In the presence of the CP- violating term f^, our analysis of CP- violating obsevables 
in Ref. ^ can be easily applied here. 

However, as the coupling of iiZ^ is almost purely axial-vectorial, there is approximate 
charge symmetry C, which assigns this vertex even C-parity, with the f^ term also even as 
well. Any C-odd observable will be suppressed. 

We find out that the most prominent effect of CP nonconservation resides in the elements 
(+, — ) or (— , +) of the spin-density matrix, 
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Ar(e) ^^^'^""' V— 7 / vA^ . A^ 7— ^ 4/52 sin^ e + 7-4 ■ 

(14) 

This particular location in the density matrix produces the azimuthal dependence in the 
form of sin 20. If we integrate ip and over quadrants, we expect that CP nonconservation 
appears in the folded asymmetry, ^"(B), which is 

[dN{i, 9, 1+III) + dN{i, TT - e, I+III)] - [dN{i, 9, II+IV) + dN{i, tt - 6, II+IV)] 

[dN{e, 9, 1+II+III+IV) + dN{£, vr - 9, 1+II+III+IV)] ' ^ ^ 

Here the range of the azimuthal angle has been divided into four usual quadrants I, II, HI 
and IV. It turns out that this observable A" is C-even and thus it is not subjected to the 
suppression from the approximate C symmetry. 
We can show that 

^"(9) = --(lmp(9)+_-Imp(7r-9)_+) . (16) 

In Fig. 1, we show the CP-odd asymmetry in the density matrix versus the scattering angle 
9 per unit of small Re // at various energies, i/i = 200, 250, and 300 GeV. Observation 
of this asymmetry is a genuine signal CP violation, as it is not faked by the final state 
interaction. 

It is interesting to note that we do not need to know the charge of i as the events are 
collected over quadrants I+III or II+IV. We can use this fact to apply our formula even to 
the larger sample of jet events from the Z° pair without tagging the charges of the primary 
partons. Our formalism can be easily translated for the process qq -^ Z^Z^ in the hadron 
collider. 



IV. TWO HIGGS MODEL 

Cubic couplings among neutral gauge bosons do not appear at the tree level in the 
standard model gauge group of SU(2)2,xU(l)y. But they can be induced at the loop level. 
In the minimal standard model with just one Higgs doublet, such amplitudes do not have 
any CP violation even at the one-loop level, as will be clear from our analysis below. We 
therefore perform the calculation of CP violating effects in these trilinear couplings when 
there are two Higgs doublets [|1^ present in the model, which is a popular model in its 



own right. Among the possibihties which open up with the two doublets are: spontaneous 



CP violation |]TT|, incorporation of the Peccei-Quinn symmetry [Q to solve the strong CP 
problem, and incorporation of supersymmetry. 

At the one-loop level, cubic coupling obviously comes from triangle diagrams. If the 
internal lines are fermions, no CP violating effect is generated at the one-loop level, because 
the Z or photon couplings with fermions are flavor diagonal and CP conserving. There 
are also triangle diagrams with internal W lines. In the Feynman-t'Hooft gauge, it can 
be shown that they do not contribute to the form factors as shown in Eq. (|I|). Thus, for 
our purpose, we need to calculate only the diagrams involving Higgs bosons in the loop. 
Obviously, such diagrams can never involve the antisymmetric e-symbol, so one can only 
obtain a non-zero //. This term has been shown to be non-zero for WW Z coupling at the 



one-loop level for the model at hand |I^. We want to extend their calculation for the case 
of V*ZZ couplings, where V* can be either an off-shell Z-boson or photon, and the other 
two Z-bosons are assumed to be on-shell. 

To set up the notation, we call the two Higgs multiplets to be ipi and ip2- Usually, 
they are assumed to have special transformation properties with respect to some discrete 
symmetries in order to avoid flavor changing neutral currents at the tree level. We assume 
that such discrete symmetries are not imposed on the soft terms in the Higgs potential, 
otherwise CP violation would be eliminated in the Higgs sector of the model. Without any 
loss of generality, we can take the vacuum expectation values (VEVs) of ipi and ip2 to be 
Vi exp(i'(9) and V2- One can then deflne a linear combination yj of the two multiplets which 



has a VEV v = Jvf + V2, and the orthogonal one, ip', has a vanishing VEV. The components 
of these doublets can then be written as 

The fields shown here are complex combinations of the fields in the ifi-if2 basis. The compo- 
nents w^ and z are eaten up by the gauge bosons and disappear from the physical spectrum. 
There are four physical spinless bosons in the model. One of them is the complex field H~^. 
The other three are, in general, superpositions of the fields 0i, 02, and 03. We define the 
eigenstates by Ha, where 

3 

(pa=T.OaAHA, (18) 

A=l 

O being an orthogonal mixing matrix. 

The coupling of these neutral Higgs bosons with the Z-boson looks very simple in the 

0-basis: 



Vertex 



Mp) -^ ^(?) 
02 (p) -^ 03(g) 



Feynman rule 



27i^(P + ?)M (19) 



Using Eq. (0), it is trivial to rewrite these Feynman rules in terms of the mass eigenstates 
of neutral Higgs bosons: 
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Vertex 



Ha{p) -^ z{q) 
Ha{p) ^ Hsiq) 



Feynman rule 
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2cose^O,A{p + q), (20) 

2^^(02a03b - 02bO,a){p + g)^ 



Using the orthogonality of the mixing matrix O, we can write 

02^03^ — 02503^1 = 2^ ^ABcOic , (21) 

c 

which simplifies the form of the Z-coupling with two physical Higgs bosons. Notice that the 
Z-coupling between two physical Higgs bosons is necessarily fiavor- changing, which opens 
up the possibility for CP violation at one-loop level. For the reason that the photon field 
preserves fiavors at the tree level, there is no fj form factor at the one-loop calculation in 
the Two Higgs Model. 

These cubic couplings appear in the triangle diagrams shown in Fig. 2. Notice that, in 
the figure, the Higgs boson lines have been denoted with subscripts z, j, fc, which run from 
to 3, where Hq is identified with the unphysical Higgs z which appears as intermediate lines 
since we adopt the Feynman-t'Hooft gauge. A straightforward calculation now shows that 
the form-factor // from these diagrams can be written in the form 

Here, \ijk is a factor coming from vertices which will be discussed below, and the loop 
integral /(Mj, M^, Mk) is equal to: 



2! y y (x y) In ^^Mi-wil-w)Ml- xvP^ - zO+ '^'^'^^ ' ^^^^ 



where the positive Feynman parameters x and y are restricted within the integration domain 
x + y < I and also w = 1 — x — y. A is a cut-off which disappears in the expression for ff , as 
we will show below. When one of the particles denoted by i, j or k is the unphysical Higgs 
boson, the corresponding mass should be interpreted to be Mz, because the propagator of 
the unphysical Higgs boson has a pole for this value of mass in the gauge we use. For future 
purposes, notice that 

/(M„ Mj, Mk) = - I{M,, M„ Mk) , (24) 



which follows from the definition in Eq. 

Let us now discuss the factor Xijk- First, consider the case when all the Higgs bosons in 
the loop are physical ones. Due to the antisymmetry of the coupling of HaHbZ^ from Eq. 
(pop, all the Higgs bosons in the loop must be different. If, following the direction of the 
momentum arrow in Fig. 2, we encounter the mass eigenstates iJi, H2 and H'^ in that order, 
it is easy to see that the factor coming from the vertices is 

Ai23 = OnOi20i3 = A. (25) 



Obviously, there are three such diagrams, and their total contribution is 

A {/(Ml, M2, M3) + /(M2, M3, Ml) + J(M3, Ml, M2)} . (26) 

On the other hand, if we encounter the eigenstates in the reverse order, we obtain a factor 
—A from the vertices. However, this term will be multiplied by 

{/(M2, Ml, M3) + /(M3, M2, Ml) + /(Ml, M3, M2)} . 

By virtue of Eq. (0), the product of the two is the same as the contribution of Eq. (p6|). 

Next we consider diagrams where one of the internal lines is the unphysical neutral Higgs 
boson z. Note that since there is no coupling of the Z-boson with two unphysical Higgs 
bosons, at most one internal line can be the unphysical Higgs boson. In this case, one can 
derive that 

A12O = -^230 = -^310 = ~^ , (27) 

and the same value for any even permutation of subscripts, but opposite sign for an odd 
permutation. Therefore, the last factor of summation in Eq. (|22D becomes 



^ A,,,/(M„ Mj, Mk) = 2A{ + /(Ml, M2, M3) + /(M2, M3, Ml) + /(M3, Mi, M2) 

- /(Ml, M2, Mz) - I{M2, M3, Mz) - I{Ms, Ml, Mz) 

- I{Mz, Ml, M2) - I{Mz, M2, M3) - I{Mz, M3, Ml) 

+ I{Mz, Ml, M3) + /(Mz, M2, Ml) + I{Mz, M3, M2)} . (28) 



One can see that the cutoff A dependence is cancelled by pairs in Eq. (pS]). We also note 
that ff remains finite when P^ = M| as noted in Ref. 0. 

Fig. 3 shows the extremely tiny size (~ 10~^) of ff for typical choices of parameters. We 
only use this Two Higgs Model as an illustration how CP violation occurs even in a purely 
bosonic sector. 



V. CONCLUSION 

At LEP II, the Z^Z^ production cross-section is about 1 pb (See Fig. 4) for y/s = 200 GeV 
which can provide about 500 Z^Z'^ pairs per year for the design luminosity of 5 ■ 10^^ cm~^ 
s~^. As we have shown in the paper, it is possible to test CP symmetry in purely charged 
leptonic, purely hadronic or mixed channels of the two Z^ boson decays. We may require 
that at least one of the Z° decays into the charged leptons in order to avoid backgrounds from 
the W^W~ production. The branching ratio of a single Z^ decaying into all charged leptonic 
channels (e+e~ + /i"*"//" + t^t~) is about 10%. While the event statistics probably will not 
be large enough to test some of the popular alternative gauge models of CP violation, it is 
nevertheless sufficient to provide nontrivial constraints on the CP-odd form factors in the 
three gauge boson couplings. 
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FIGURE CAPTIONS 

Fig. 1 The CP-odd asymmetry in the density matrix versus the scattering angle per unit 
of Re // at various energies, y/s = 200, 250, and 300 GeV. 

Fig. 2 Triangle diagrams with internal scalar lines which give rise to the Z*ZZ coupling. 

Fig. 3 The size of ff /OuOuOis versus the lightest Higgs mass at ^/s = 200 GeV, for the 
case M2=150 GeV, and M3 = 250 GeV. The real and the imaginary parts are given 
by the solid and the dashed lines respectively. 

Fig. 4 Differential cross-section da/dcosQ for e^e+ -^ Z'^Z^ at various energies ^/s= 200 
(solid), 250 (dashed), and 300 GeV (dased-dotted), predicted by the Standard Model. 
The horizontal lines indicate the level of the corresponding total cross-sections. 
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